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Abstract 

Measurements of resonant tunneling through a localized impurity state are 
used to probe fluctuations in the local density of states of heavily doped 
GaAs. The measured differential conductance is analyzed in terms of corre- 
lation functions with respect to voltage. A qualitative picture based on the 
scaling theory of Thouless is developed to relate the observed fluctuations to 
the statistics of single particle wavefunctions. In a quantitative theory correla- 
tion functions are calculated. By comparing the experimental and theoretical 
correlation functions the effective dimensionality of the emitter is analyzed 
and the dependence of the inelastic lifetime on energy is extracted. 

PACS numbers: 73.23.-b, 72.20.My, 85.30.Mn 
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I. INTRODUCTION 



The observation of impurity-assisted tunnehng in vertical transport experiments on 
double-barrier semiconductor heterostructuresS^ led to the possibility of using the reso- 
nant impurity level as a local probe of electronic states of electrodes prepared from heavily 
doped degenerate semiconductors. A number of experiments in strongly asymmetric double- 
barrier structures have measured directly the local density of states (LDOS) of an electrode 
as a function of excitation energy E from the Fermi level,iHl^ those with the highest spec- 
tral resolution reporting features including Zeeman splitting of single-particle levels in a 
disordered emitter.! 

The idea of such experiments is illustrated by the sketch in Fig. 0. Electrons tunnel from 
a heavily-doped disordered emitter through the energetically lowest level of the quantum 
well sandwiched between the double barriers. This energetically lowest level of the quantum 
well serves as the spectrometer S. At zero bias, the energy of this impurity level, Es, 
does not coincide with the chemical potential fi in the emitter. It comes to resonance 
only after the bias voltage reaches a threshold value Vs{Es). Typical current-voltage I(V) 
characteristics of such a device can be divided into three intervals:i^'0 one interval below 
the threshold, where / ~ 0; the threshold regime V ~ Vs{Es), where I{V) undergoes a 
jump when the resonant level crosses the Fermi level /i in the emitter, and the interval 
of a plateau, Vs{Es) < V < Vi{Ei), where the current remains nearly constant. The 
latter interval lasts until the next impurity level Ei is lowered enough to contribute to the 
transport and it is ideal for studying the image of the LDOS in the emitting reservoirjllli 
since any further variation of the current as a function of bias voltage, I{V) is dominated 
by the energy dependence of the tunneling density of states in the emitter, I{V) oc ^{E). 
A convenient way to look at such I(V) characteristics is to plot the differential conductance 
G(y) = dl/dV oc du/dE, in which the image of variation in LDOS is more pronounced. 

In a disordered medium, the energy dependence of the LDOS studied at a certain point 
of a sample reveals an irregular fine structure,^! ^{E) = 1^0 + ^^^i^), which arises from 
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quantum interference of elastically scattered quasi-particles diffusing coherently within a 
length scale related to their lifetime at a particular energy. Since the Aharonov-Bohm phase 
accumulated by a diffusing particle in a magnetic field changes the interference pattern, such 
a fine structure, 5v{E) depends randomly on a magnetic field B. In tunneling experiments, 
the interfering quasiparticle is, in fact, a 'hole' in the Fermi sea left behind by the tunneling 
of an electron out of the emitter with E < fi. Being in a non-equilibrium state, such a hole 
'floats up' towards the Fermi level, due to inelastic collisions between electrons, so that it 
can be characterized by a finite lifetime equivalent to a broadening of emitter states.0 The 
broadening of emitter states washes out the finest features in the LDOS fluctuations, and, 
therefore, it strongly affects the amplitude and correlation parameters of fluctuations of the 
differential conduct ancelll of a given resonant tunneling device, SG{V) = G{V) — (G). A 
particularly convenient situation to study fluctuations is realized in devices where the mean 
value of the density of states in the emitter and also the transmission through the barriers 
varies much slower than fine fluctuations in LDOS, so that within the narrow energy interval 
below the Fermi energy of the emitter (G) is negligible and 6G{V) ~ G{V). 

Recently, we reporte an experiment where the speeding up of the quasi-particle relax- 
ation upon the increase of the energetic distance to the Fermi level (equal to the excitation 
energy of the Fermi sea holes) was observed via the decline in the variance of differential 
conductance fluctuations at higher bias voltages, ((^G)^). In the present publication, we 
study the correlation function, K(y), of a random differential conductance pattern worked 
out for different bias voltage intervals. To make this analysis sound, the differential conduc- 
tance of a resonant tunneling structure has been measured for a dense grid of magnetic field 
values, which has largely increased the statistical ensemble of data used in the evaluation of 
the correlation function K{AV) of a random pattern 6G{V,B), 

= , (1) 

and has allowed us to compare details of its shape to the results of a theoretical analysis. 
The latter effort has enabled us to notice some geometrical features of the structure used 



in this experiment (produced in a particular growth process), which would be difficult to 
detect otherwise. 

The material presented in this paper is organized as follows. In Section II, the exper- 
imental set-up, parameters and design of the structures which we used, and the raw data 
of G(V) characteristics are discussed. In order to characterize the spectral resolution of the 
spectrometer, we analyze in Section II the form of {G(y))B averaged over many runs taken 
at different values of applied magnetic field. A detailed quantitative analysis of fluctuations 
and their correlation functions is presented in Section III, in comparison to the results of a 
theory presented in two Appendices. The end of Section III is devoted to the discussion of 
energy dependence of the quasi-particle relaxation rate extracted from this analysis, from 
the point of view of Aronov-Altshuler theory of electron-electron interaction in disordered 
metalfl0. Appendix A completes the text with a qualitative estimation of the variance of 
the differential conductance fluctuations based upon the theory of statistical and correlation 
properties of chaotic wavefunctions in disordered media using an approach similar to Thou- 
less's scaling theory.lil The quantitative analysis of the variance and correlation properties 
of a pattern of 6G{V), including the dimensional crossover, is presented in Appendix B. 



II. SAMPLES AND EXPERIMENTAL RESULTS 

The experiment was performed using an asymmetric double-barrier heterostructure 
which was grown by molecular-beam epitaxy on n+-type GaAs substrate. Directly on top of 
the substrate the layer sequence for the resonant tunneling diode was grown as illustrated in 
Fig. H. The growth started with a 300 nm thick GaAs layer doped with Si to 4.0 x 10^'^ cm~^. 
This emitter layer is followed by a very thin spacer layer of 7 nm undoped GaAs. The actual 
resonant tunneling structure consists of a 10 nm wide GaAs quantum well sandwiched be- 
tween two Alo.3Gao.7As barriers of 5 and 8 nm width (top and bottom barrier). The collector 
of the structure is formed by a second spacer layer of 7 nm undoped GaAs and a 300 nm 
thick layer of GaAs doped with Si to 4.0 x 10^^ cm~^. 
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The barrier structure is highly asymmetric, the transparency of the thick emitter barrier 
is much lower than that of the collector barrier, which means that the value of the tunneling 
current is dominated by the low transmission of the emitter barrier. Due to the thin spacer 
layer the nominally undoped quantum well contains a small number of residual impurities. 
The energetically lowest impurity state will be used as a local spectrometer of the emitter 
states. 

In order to limit the number of residual impurities in the quantum well, pillars with small 
areas were fabricated from this heterostructure.0 By employing electron-beam lithography, 
evaporation and lift-off, AuGe/Ni layers were deposited on the top of the wafer. This 
metalization served both as an ohmic contact and as an etch mask for the following reactive 
ion etching (RIE) step. A AuGe/Ni coating was also evaporated onto the substrate side 
of the wafer to form the back ohmic contact. Free standing pillars with diameters in the 
/im and sub-/im range and a typical height of several hundred nm were etched using RIE. 
Then, large-area Cr/Ag/Au bond pads could be prepared on top of the pillars by planarizing 
the pillars with an insulating polyimide layer. The tunneling current was measured with 
a dc technique in a dilution refrigerator at 20 mK base temperature. For our analysis the 
differential conductance G{V) was numerically calculated from the measured current values. 
A typical G{V) trace is shown in Fig. ^ At zero bias, S lies above the Fermi level in the 
emitter and is not available for resonant transport, resulting in G = 0. At Vs = 9.8 mV, the 
spectrometer crosses the Fermi level and the current jumps abruptly from zero to a finite 
value, resulting in a sharp peak in the derivative G oc du/de. For larger bias voltages a 
reproducible oscillatory fine structure can be seen, which we attribute as the result of LDOS 
fluctuations in the contact regions. This fine structure is formed by electrons which tunnel 
from below the emitter Fermi level through the lowest discrete state S in the quantum well. 
Since the emitter barrier is much stronger than the collector one, the value of the current 
step is mainly determined by the tunneling rate F;/^ through the thick barrier on the emitter 
side. Due to this large barrier asymmetry the G(V) curve at voltages V > 9.8 mV represents 
the energy dependence only of the LDOS dv/de in the emitter contact. So the fine structure 
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represents an image of emitter contact LDOS fluctuations scanned by the impurity-related 
level in the quantum well.l 

The quantum interference interpretation of the observed fine structure is supported by 
the observed effect of an applied magnetic field. The oscillatory form of G{V) randomly 
changes upon variation of a magnetic field, at the scale AB < 30 mT. Figure ^ shows a 
grey-scale image of the differential conductance measured as a function of both bias voltage 
and a magnetic field within the interval of fields — lT<i?<lT, where the Landau 
quantization of states in the emitter is completely suppressed by disorder. This diagram 
is symmetric with respect to magnetic field inversion, as it should be for a two-terminal 
measurement. The use of magnetic field enables us to get a sound amount of data for the 
following statistical analysis of fiuctuations. 

The onset of resonant tunneling through the lowest lying impurity state S appears in 
Fig. ^ as a black line at a voltage of 9.8 mV (in parallel to the B-axis). The second black 
line at a bias voltage of 14.6 mV appears when the next, higher lying impurity state crosses 
the Fermi energy in the emitter. For voltages ranging from 9.8 mV up to 14.6 mV, the 
measured tunneling current results only from tunneling through the lowest lying impurity 
state S. This state is used as a local spectrometer to scan the LDOS below the Fermi 
level in the emitter. Voltage and energy scales are related via E = ae{V — Vs) where 
the prefactor a = 0.5 accounts for the fact that only part of the voltage drops between 
emitter and spectrometer Therefore, the plot in Fig. ^ covers an energy range of quasi- 
particle excitations of about < E < 2.4 meV, which is indicated by an alternative scale 
for the horizontal axis on the upper side of this figure. The amplitude of G{V) fiuctuations 
decreases (fine structure washes out) with increasing bias voltage in the range 9.8 mV < V < 
14.6 mV (interval between two peaks). At the same time, the characteristic voltage scale 
dominating the fine structure increases which is interpreted below to be the result of the 
inelastic broadening of quasi-particle states in the emitter. Note that although oscillations 
at larger energy scales are also present in I{V) oc v{E), their contribution to G{y) is 
suppressed due to the differentiation. For a broad spectrometer both the amplitude and 
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correlation voltage of fluctuations would be the same over the entire range of Vs' < < Vi. 
For a narrow spectrometer, as studied in the present work, inelastic broadening of states in 
the bulk exceeds the spectrometer width upon increasing the excitation energy of a quasi-hole 
left in the emitter. Then, this inelastic broadening affects the parameter of the fluctuation 
pattern. Note that the observed fluctuations become sharp and large again after the second 
impurity level begins to contribute to the current aXV> 14.6 mV. This is because tunneling 
through the second impurity state involves states close to the Fermi energy which have 
negligible inelastic broadening. In the following, we shall focus on the tunneling through the 
lowest lying impurity state, i.e. on the interval of bias voltages smaller than V < 14.6 mV. 

In Fig. ^the differential conductance is shown after averaging the raw data G{V, B) over 
the interval of magnetic field specified above. This averaging increases the contrast between 
the main peak corresponding to the spectrometer S crossing the emitter Fermi level and 
(Giy)) g at larger bias voltages, where a random contribution from LDOS fluctuations is 
strongly suppressed. The fluctuations are suppressed by a statistical weight of -\/iV, where 
is the number of uncorrelated G(V) traces taken at various magnetic fields. The plot in 
Fig. ^ can be used to extract the nominal spectrometer width, F. The emitter barrier in 
the device we study is thicker than the collector one, so that the broadening of the resonant 
level is dominated by electron escape from it to the collector, F = F^, + F; ^ F^, whereas the 
value of the current step is mainly determined by the tunneling rate Ti/h through the thick 
barrier on the emitter side. The averaged {GiV)) ^ characteristics at the threshold can be 
parameterizecl0^0 by the height of the conductance peak at the threshold voltage Vs and 
by its width Vr at the half-maximum, which is given by,ii 

Vr^T/{ea). (2) 

Below, we use F = 36 /ieV taken directly from Fig. ^ 
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III. STATISTICAL CHARACTERISTICS OF THE DIFFERENTIAL 
CONDUCTANCE FLUCTUATIONS PATTERN. 



In this section, we analyze correlation properties of the measured differential conductance 
pattern, aiming to extract from this the value and the form of the energy dependence of the 
decay rate for quasi-particles. The relevance of the correlation function of the fluctuations 
pattern for such an analysis arises because the auto-correlation function of fluctuations re- 
flects the typical scale of their energy dependence (which is equivalent information to that 
in the power spectrum of frequencies of oscillations). To start with, the pattern of 5G{V) 
is random and it is related to the derivative of the LDOS with respect to energy, where the 
contribution from features at the flnest energy scale is enhanced by differentiation. There- 
fore, the correlation function K{AV) = {SG{V+AV, B)5G{V, B))/{SG^) carries information 
about the flnest resolution of quantum states in the emitter. On the one hand, due to the 
flnite spectrometer width, F, flne structure in the LDOS at energy scales smaller than F 
is smeared out by the spectrometer, so that in the measurements reported above it cannot 
be resolved. On the other hand, the flnest energy scale of LDOS fluctuations is intrinsi- 
cally limited by inelastic broadening of quasi-particle states in the emitter, h^y. As a result, 
the typical value of bias voltage, at which the differential conductance varies randomly, is 
determined by the sum of the above two, 

V,= -[r + hy] = E,/a, (3) 
a 

where the spectrometer width F is the same for the entire interval of energies of the quasi 
hole in the emitter (left behind by the tunneling process) that we are able to study using 
one impurity state, whereas the inelastic broadening, h'y{E), is dependent on the excitation 
energy and varies across the studied bias voltage interval. The same combination of energetic 
parameters also determines the variance of the differential conductance fluctuations, {SG"^), 
which will be discussed in subsection B. 
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A. Correlation function of fluctuations. 



The experimental determination of the correlation function, K, consists of the evaluation 
of the variance {5G'^)b and, then, the autocorrelation function of the measured differential 
conductance fluctuations pattern, 

K{^V) = {5G{y + /\V,B)5G{y,B))B/{5G^)B. 

by means of averaging over different magnetic field points within the interval 0< B <1T. 
Then, the obtained correlation function is additionally averaged over a narrow interval of 
bias voltage, not more than 2-3 times broader than the width of the autocorrelation function 
determined after the first step. This procedure allows improved statistics and it slightly 
reduces variations in the form of the correlation function. Note that the finite amount 
of data used in this analysis still leaves space for statistical errors, so that the evaluated 
correlation function may be treated seriously only within an interval equal to 3 times its 
width at the half-maximum. 

The typical result we get for such a correlation function is shown in Fig. ^ for two values 
of bias voltage: one at the beginning of the studied interval, at = 10.2 mV, the other - at 
its end, at ^ = 13.8 mV. These correlation functions have a very different width, which we 
attribute to an increase of inelastic broadening of states of quasi holes in the emitter upon the 
increase of their excitation energy, such that it becomes even larger than the spectrometer 
width, r. Therefore, the comparison of correlation parameters of K{/SV) can be used for 
determining directly the value of the inelastic relaxation rate of quasiparticles in the emitter 
function of their excitation energy. 

To obtain an absolute value of the inelastic broadening from such a comparison, one has to 
make a certain fit and, therefore, to use a certain form of the correlation function K{^V). 
Theoretical analysis of the correlation function of differential conductance fluctuations in 



Ref. 12 has shown that its form, and, therefore, the value of the correlation voltage extracted 



from the fit depend on the effective dimensionality of the diffusive emitter, that is, on its 
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geometry. In particular, for a quasi-OD emitter (diffusive pillar) and a quasi-2D film we have 
calculated 



Ko{AV) 
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For a quasi-lD wire, and 3D bulk, these are, respectively, 



Y = ^l + {AV/V^f. (5) 

All these correlation functions were obtained in the unitary symmetry class limit for fluctu- 
ations. 

In Fig. ^ all four of them are compared to the experimentally determined correlation 
function for the smallest bias voltage interval, i.e., for V = 9.8 mV. Theoretical curves shown 
in this plot for various models of an emitter can be characterized by the depth of a negative 
anti-correlation overshoot in K, which is the most pronounced in the quasi-OD case. For 
each theoretical curve, the fit to the data is made using a single parameter, Vc, and the best 
agreement between the theory and experimental data is achieved for the quasi-OD model of 
the emitter. 

The suggestion that the emitting electrode in the studied structure has the form of a box, 
rather than the form of a wire which would be a natural assumption based upon the shape 
of the lithographically processed material in Fig. 0, needs an explanation. The point is that 
the emitter side of this device has been produced by overgrowing heavily doped GaAs:Si 
substrate (10^^ cm~'^ of Si) with a 300 nm buffer layer of GaAs:Si 4 x 10^'' cm~^. It is known 
that the interface between the substrate and the first grown layer is not as perfect as the 
interfaces produced during the MBE growth process. It is expected that at the interface 
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between the substrate and the first layer a higher density of background impurities are 
incorporated and also that the dislocation density will be higher than in the rest of the 
structure. The predominant background impurities will be carbon impurities, which act as 
acceptors in GaAs and compensate the Si-donor doping. Due to this compensation this 
interface could be poorly conducting. Although poor conduction through this interface does 
not affect the observable resistance value of the device and no other measurement performed 
on structures from the same series had enough sensitivity to indicate its presence, the LDOS 
fluctuations measurements appear to be sensitive enough to illuminate its existence. 

For each given geometrical shape of the emitter (for this sample, a L-thick disk with a 
radius R), the effective dimensionality reflected by the shape of the correlation function in 
Eq. (H) also depends on the ratio between the diffusion length. 



a disordered system that is effectively tested by a coherently diffusing particle within the 
time scale taken before it either escapes from the emitting electrode to the collector via 
the resonant impurity S or relaxes inelastically into states at different energies. When the 
latter length scale is the largest, ^ L, R, the correlation function of fluctuations has the 
quasi-OD form. When R > L^. > L, the flnite radius of a pillar would not matter, and the 
correlation function would have the quasi-2D form. Similarly, L > Lc > R would correspond 
to the quasi-lD result in Eq. (|). Finally, one would have to treat the regime of R,L ^ L^. 
as the three-dimensional one. 

The value of the correlation voltage extracted from the flt of the experimental data in 
Fig. 1^ using the quasi-OD model, = 80 yuV is very close to the width of the main resonance 
peak in Fig. ^ determined by the intrinsic spectrometer width, Vr = 72 /iV. Comparison of Vp 
with other values of Vc obtained from flts of experimental K{y ) using other dimensionality 
assumptions \Vc = 65 /iV, 51 /iV and 33 fiV for the QID, Q2D and 3D models, respectively] 
gives an additional argument in favor of the view that we deal here with a quasi-OD emitter. 




(6) 



and geometrical sizes, L and R. The diffusion length in Eq. (^) characterizes the volume of 
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At the same time, the relevant diffusion length Lc calculated as Lc = yDh/Ec ~ ^JDh/V 
is longer than both the sample diameter and the width of the buffer layer, which would be 
consistent with an assumption that the interface is an obstacle for electron escape to the 
substrate. 

Since the length scale Lc in Eq. (P) shortens, due to faster inelastic relaxation as the 
quasi-particle excitation energy increases, the effective dimensionality of the system may vary 
across the bias voltage interval we study. Since the sample used here has R> L, a crossover 
may take place between the quasi-OD and quasi-2D form of the correlation function that 
should be used for fitting the data in the broader voltage interval. One can find indications 
of such a crossover in the series of correlation functions shown in Fig. ^. 

Traces of crossover behavior in Fig. ^ require one to make a detailed theoretical analysis 
of the intermediate regime Lc ~ R, since our final goal is to obtain quantitative informa- 
tion about the quasi-particle lifetime, as a function of quasiparticle energy in the entire 
energy interval assessed in the reported measurement. Details of a calculation of correlation 
functions in the crossover regime are presented in Appendix B. Here, we only describe the 
results, in a graphic form. Fig. ^ shows the change in the shape of the correlation function 
of differential conductance fluctuations expected for a spectrometer placed in the center of 
the bottom surface of a round disk, for various values of the ratio Lc/R, but for the same 
nominal Vc- This plot shows that the crossover between the quasi-OD (dotted bottom line) 
and quasi-2D form (solid line) can be split into two steps. First, the negative deep in K{V) 
aXV^Vc is reduced (anti-correlations become weaker), which happens without a noticeable 
change in the width of the correlation function at the half maximum (in units of V/Vc). The 
following evolution of the form consists of the broadening of the main part of the correlation 
function. This two-step evolution suggests that the fit to the central peak of the experimen- 
tally determined correlation function using the quasi-OD formula is a consistent procedure 
applicable even across some part of the crossover regime. The need for such a simplified 
procedure in the following analysis has another reason. When the crossover takes place, the 
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exact form of K{V) becomes dependent on the position of the spectrometer on the surface, 
that is, its distance to the disk perimeter. This effect is illustrated in the inset to Fig. | using 
several plots of K{AV) calculated for different off-center positions of the resonant impurity. 
Plots in Fig. ^ also show that values of the sum of the physical parameters, F + ^7 obtained 
following such a procedure may be overestimated when the crossover to the quasi-2D limit 
is more developed. 

B. Analysis of the variance of differential conductance fluctuations 

Quantitative information about the energy dependence of inelastic quasi-particle relax- 
ation can also be extracted from the bias voltage dependence of the variance of differential 
conductance fluctuations. Such a dependence for the sample described in this paper is shown 
in Fig. It is evaluated on the basis of the pattern of raw data in Fig.^ after subtracting 
from the data the average conductance, G{V), shown in Fig.^, then, averaging the difference 
over the magnetic field interval 0< i? <1T 

varBG=/ —{G{V,B)-{G{V))sf, (7) 

JOT 1 i 

and, then, by smoothing it over the bias voltage interval of 3 times Vc determined for the 
corresponding bias voltage range in the previous section. The result is presented in the form 
normalized by the height of the main conductance peak. Or, in order to exclude from this 
analysis the parameters of tunneling barriers, and the bias voltage value is converted here 
into the excitation energy of a quasi-particle {E is the energy of the Fermi sea hole evaluated 
with respect to the Fermi level). Because of the above mentioned smoothing procedure, we 
cannot start the plot in Fig. |10| from exactly E = 0. 

The decrease of the amphtude of differential conductance fluctuations upon the increase 
of excitation energy of quasi-particles is attributed to a faster inelastic relaxation of the lat- 
ter, which can be used to study the dependence 'y{E). Similarly to the correlation function, 
the exact form of such a dependence varies if one makes different assumptions about the 
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effective dimensionality, d: 
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Using the measured spectrometer width F and the known sample dimensions R and L, 
these equations enable us to obtain theoretical estimates of the amplitude of the variance 
for a given effective dimensionality. A comparison with the low E part of the measured 
variance data plotted in Fig. |10|, where we expect 7(-E) ~ 0, has the best agreement with 
the QOD theory. In these estimations we used the value of the mean free path, / ^ 70 nm, 
assigned to the nominal doping level of the buffer layer. This value of / is confirmed by the 
tendency of the variance {6G{By) to follow a [1 + (co'cT)^]-dependence@'i at classically high 
magnetic fields uJcT ~ 1 as shown in the inset of Fig. p!0| . 

The energy dependence of the parameter = T + h'y can be extracted from Fig. ^ 
using the formulae in Eq. (|^) and it is plotted in Fig. |TT] for the four different effective 
dimensionalities (upper, solid lines). Also plotted in Fig. |Tl| is the energy dependence of the 
parameter Ec = aVc obtained from the analysis of the correlation function (lower, dashed 
lines). The values of Ec obtained along two different roots have to coincide for an appropriate 
dimensionality assumption, and they agree only when analysis is based upon the quasi-OD 
emitter model. 



C. Quasi-particle inelastic relaxation rate in a disordered conductor. 



On the basis of the material presented above, we conclude that the use of the quasi-OD 
assumption for the analysis of fluctuations is fully justified and can be exploited for analyzing 
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the energy dependence of the inelastic relaxation rate of quasi-particles, '~f{E). The latter 
can be obtained from the data shown in Fig. |Tl] by subtracting the original spectrometer 
width. The resulting relaxation rate dependence on the excitation energy is shown in Fig. |T^. 
This plot contains two sets of data taken from the analysis of correlation functions and the 
variance, and the comparison to the rate values calculated using Altshuler-Aronov theory. 
The discrepancy between data worked out in two different ways indicate the arrow bars one 
would have to assign to the presented analysis. 

The theoretical curve shown in Fig. |12] is a fit to the relaxation rate as derived by Sivan, 
Imry and AronovJ^^ using Ep = 30 meV for the emitter buffer doped to 4.0 x 10^'' cm~^ 
with Si, 

^ ~ 167r r3/2E2 ■ 

The mean free path obtained from this fit is / = 93 nm, which is close to mean free path 
expected for this nominal doping (between I = 50 nm and / = 100 nm) and also close to 
the value extracted from the analysis of the increase of the variance of fluctuations with 
magnetic field (/ = 70 nm). The use of the three-dimensional expression for the relaxation 
rate in Eq. (^, in contrast to the quasi-OD model used to describe fluctuations, is justi- 
fied by the following reason. As discussed at the end of Appendix A, the relaxation of a 
quasi-particle with energy E is dominated by electron-electron (e-e) collisions with energy 
transfer comparable to E, and such a rate is determined by correlations between chaotic 
wavefunctions with a typical energy separation e ^ E. The latter are formed at the length 
scale ~ ^hD/e, which has to be compared to the system size: the pillar radius, R, and 
the width, L. In particular, if < L < R, the e-e interaction can be treated as in the 
three-dimensional bulk of a disordered conductor. This condition can be expressed more 
rigorously as E > tt'^HD/L'^, -k^HD / R? which states that the quasi-particle excitation energy 
has to be larger than the Thouless energy related to diffusive motion across the pillar. Since 
the extracted values 7(-E') sufficiently exceed experimental uncertainty only for quasi-particle 
excitation energies E > 0.5 meV [which has to be compared to n'^hD/L'^ ^ 0.4 meV], their 

15 



quantitative comparison to the calculation of 'y{E) in the 3D limit seems to be consistent. 



At the same time, the entire interval of energies analyzed in Fig. ^ belong to a clearly 
diffusive regime, E < hr ^ A meV. Note that all this analysis is extended only over the low 
magnetic field range, where the Landau quantization of emitter states does not play any 
role. 



IV. SUMMARY 

We study resonant tunneling through a discrete localized level in a GaAs/AlGaAs double- 
barrier heterostructure. The differential conductance exhibits a temperature-insensitive fine 
structure which is attributed to fluctuations in the local density of states in the doped GaAs 
emitter. The observed fine structure is analyzed in terms of variance of the fluctuations in the 
differential conductance and in terms of correlation functions with respect to voltage. From 
analyzing the shape of the correlation function we conclude that the effective dimensionality 
of the emitter is zero-dimensional caused by the disordered interface between the GaAs 
substrate and the doped buffer layer. In this experiment the electrons tunnel from below 
the Fermi energy in the heavily doped emitter contact through the discrete localized level 
leaving behind a quasi-hole in the emitter. By quantitatively analyzing the width of the 
measured correlation functions and the measured variance we are able to extract the energy 
dependence of the inelastic quasi-hole relaxation. 
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APPENDIX A 



This appendix presents a quahtative method of estimating the variance of differential 
conductance fluctuations and the energy dependence of quasi-particle relaxation. It is con- 
structed using a scaling picture similar to that of Thoules^l by considering what happens 
to the states of single electrons in a box when the electrons are able to diffuse into other, 
similar boxes. For clarity we begin by considering three dimensions d = 3 although this 
is not necessary for the following arguments to hold. In a classical picture of diffusion, a 
diffusive path can be viewed as a series of straight line segments of typical length equal 
to the elastic mean free path /, where / = vpT, t is the elastic time, and vp is the Fermi 
velocity. The classical diffusion coefficient is D = Vpr/d and the typical time required to 
diffuse a length ^ is td{^) = ^'^/D. 

Consider eight cubes of length ^ which are separated by barriers such that no particles 
may move between the cubes. We imagine that it is possible to diagonalize the Hamiltonians 
of the separate cubes and we denote the eigenstates as ipai{r) where a specifies the cube 
and i specifies the state. These states are called 'mother states' of the generation ^ and they 
have a mean level spacing A(0 where A(0 = l/iK"^)- The states are normalized so that 
/ \i^ai{^)\'^dr = 1. The Hamiltonian of the total system, consisting of eight cubes of size ^ 
with barriers between them, is also diagonal. 

When the barriers between the cubes are removed, particles may diffuse between them. 
As stated above, the typical time to diffuse a length ^ is T£)(^) = ^"^/D. The energy cor- 
responding to this time is called the Thouless energy E{C,) ~ hD/^^. Diffusion between 
cubes produces a finite mixing of states from different Hamiltonians so that the Hamilto- 
nian of the new system (consisting of the eight smaller cubes) is not diagonal. Instead it 
has finite elements within a distance E{^) of the main diagonal and has elements which are 
approximately equal to zero elsewhere. 

An approximation is used to diagonalize the new Hamiltonian. An area of width E(^) 
is centered on the middle of the Hamiltonian and a unitary transformation U is applied to 
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diagonalize it, neglecting the rest of the Hamiltonian. Approximate eigenstates of the new 
system are hnear combinations of a finite number of mother states of the generation i^, 

V^?n(r)^E««^i.W (10) 

ai 

where a^" are coefficients with indices ai that refer to the original cubes of scale ^ and 
index n of the new states in the cube (3 of scale 2,^. The new approximate eigenstates are 
normalized so that 



Ei«^rr-Ei«^rr = i (11) 

n ai 



where we used the property WU = 1 (only one value of j3 is considered). Correlations 
between local densities remain important because at each level in the scaling procedure 
there is only a finite basis involved in the construction of new states. 

Consider an experimental observation of the local density of states at position tq. For- 
mally the local density of states may be expressed in terms of a summation of states. When 
observing through a spectrometer of energy width F then this may be written as a sum of 
states with energy Ei within F of the spectrometer energy Es, 

u{vo;Es)^T-' E |V^L(ro)f; ^ < ^r- (12) 

\E,-Es\<T 

The approximate eigenstates ijj^^{ro) are those in a system of size ^ where ^ < and Lr is 



the length scale corresponding to energy F, Lr ~ yhD/T. In general, however, the system 
is larger than Lr and it is necessary to know how the above summation behaves at larger 
scales C, > Ly- Consider for example a system of scale C, = 2Lr. Applying the above scaling 
procedure we may write the approximate eigenstates at larger scales using Eq.(|iy), giving 

u{ro;Es)^T-' E Hi'^'i^o/ 
\E,-Es\<r 



r-^EEK"l' CHro) 



2 



n ai 

2 



r-'ECHro) (13) 



where in the last step we used the normahzation condition given in Eg. ([TTD . This result 
shows that the summation over the energy interval F will not vary when the spectrum is 
modified into its final form at the total system size L, but it will depend on the spectrum 
at length scale Lp. This is because once the scale ^ > Lr then the corresponding energy 
-E'(^) < r and information about correlations that is carried by the mother states will remain 
in the energy interval F no matter how large ^ becomes. Now we describe the application 
of the above scaling picture to the differential conductance^'! The current in the plateau 
regime, I, is determined by a sum of local densities of the wavefunctions |'?/'£;(r)|^ with 
energy, E, taken in an energy interval F around the energy Eq, / oc ~ F~^ J2 The 
number of states in a sample of volume L"^ within the energy interval F is A^(F, L) ^ i^qTL'^ 
where is the mean density of states per unit volume, per unit energy. The variance of the 
differential conductance, (5G^), is given by a typical fluctuation in the density of states 6u 
divided by the typical energy interval F. As described above, a summation over the energy 
interval F depends on the spectrum at length scale with number of states A^(F,Lr).0'^ 
Since (r) from a single state is a random variable with mainly Gaussian statistics in the 
metallic regimeil and the variance, {5G^), is given by a sum of the individual variances, we 
have 

\ I F2Vi? 
where the typical density of a single state is <^|-?/'£;(r)|^^ ~ l/-^r- 

In the above estimation Vr is the smallest voltage step which is given by the spectrometer 
width Vr ~ F/e. An additional level broadening, hq^ takes into account relaxation processes 
in the bulk of the emitter and results in a total level broadening F (1 + ^7/F). Counting 
powers of F in the estimation of (^G^) at the end of the last paragraph leads to a factor 
of {X^n^ivf'^-^ in the variance. However we should stress that the level broadening hq 
does not influence the average value of the current in the plateau regime, 

Ar(F,Lr)(|^i.(r)''' 



or the height of the main differential peak which is given by the mean current divided by 
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the width of the peak, 



NiT,Lr)(\Mr: 
Gr ~ 



|2 



TVr 

Thus the variance is parametrically reduced as compared to Gp by a factor l/A^(r,Lr) 
where N{T,Lr) ~ g{Lr) ~ z/D^/^F^^'^/^. When the dimension d refers to the effective 
dimensionahty of the system as determined by the volume over which mesoscopic fluctuations 
occur, embedded in a nominally three dimensional space, then the factor A^(F, Lr) is replaced 
by A^(F, Ly)L^~'^. In this case the variance may be written as 

1 Gl 



It is also possible to explain the energy dependence of the inelastic scattering rate in 
terms of the above scaling picture. The rate is determined by a collision between four 
particles involving transferred energy u. There are two initial particles with energies E > 
and e' < and two final particles with energies E — u > and e' + a; > 0. The inelastic rate 
may be estimated using Fermi's Golden 

0<i^<£; -uj<e'<0 ^Is/ 

where M{E, e', u) is the matrix element for the collision. For a short ranged interaction the 
matrix elements are given by a spatial integration of a product of four single particle wave 
functionJUm 

The normalized wavefunctions in a disordered system exhibit random spatial oscillations, 
each typically contributing \ip\'^ ~ After disorder and spatial averaging the product of 
four wavefunctions is roughly 



In the limit N{E{^),^) — > oo, there is an infinite basis involved in the construction of new 
states at each level of the scaling process, leading to an absence of correlations between 
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different eigenvectors. For finite N{E{^),^), however, there is a finite basis and correlations 
exist. Integration over the hypercube provides an additional factor of C,'^ so that a typical 
value of the matrix element is 

M 



Since each summation with respect to energy contributes roughly E/A{^) we find 

7(0 



iv2(£;(0,{)A(o 

Now we consider what occurs when the hypercubes are scaled up to the total system size 
L. For small energies, E < E{L) = E^h, the system is in the zero dimensional limit whereby 
the system size is always less than the length scale associated with the energy, Le, where 



Le ^ yhD/E. The above estimation holds for all scales up to the system size so that the 
inelastic rate E'^ / {g'^ (L) A{L)) where g{L) ~ N{ETh, L). For our case of interest, 

however, we sum over states up to energies greater than the Thouless energy Eth of the total 
system, E > E{L) = Eth- When the system size reaches Le, the energy scale E contains 
all the information about correlations between the states. As for the calculation of (^G^) at 
the length scale Lp, further scaling does not change the evaluation of the relaxation rate. 
The summation over the energy interval E will not vary when the spectrum is modified into 
its final form at the total system size L, but it will depend on the spectrum at length scale 
Le with number of states N{E, Le)- We replace ^ in the above estimation with L^, giving 

N\E,Le)A{Le) 
in agreement with the prediction of Altshuler and Aronov.@ 



APPENDIX B 

This appendix describes a numerical calculation of the correlation function in a disk of 
width L and radius R where L <^ R. The crossover regime between the quasi-OD {L ^ 
R -C Lc) and quasi-2D {L <^ Lc R) limits is studied. Using standard diagrammatic 
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perturbation theory techniques, it is possible to express the correlation function in terms of 
diffusion propagators in the disk.lll Since L <^ L^, the zero mode dominates the diffusion 
propagator in the direction parallel to the current flow, across the width of the disk, but it 
is necessary to sum all harmonics of diffusion perpendicular to the current flow, in the plane 
of the disk. Adopting circular cylindrical coordinates r = (p, 0, z), we consider the resonant 
impurity to be positioned at one side of the disk z = and at an arbitrary radius < p < -R 
from the center of the disk. 

The numerical procedure outlined here is necessary only in the crossover regime since the 
exact geometry of the emitter is not relevant in the limiting cases: in the quasi-2D limit a 
diffusing electron typically does not reach the boundary of the disk and it is thus possible to 
integrate over all harmonics of the diffusion propagator in the plane of the disk whereas in 
the quasi-OD limit the zero mode in the plane of the disk is not damped very effectively and 
it dominates, enabling one to neglect all higher harmonics. These approximations produce 
the analytic results given in the main text and in both cases the position of the resonant 
impurity p is irrelevant. However the position of the resonant impurity is crucial in the 
crossover regime because it is necessary to sum over many harmonics that are influenced by 
the exact geometry of the emitter. 

The correlation function of differential conductances can be obtained from the disor- 
dered averaged current-current correlation functionlli {6I{V)6I{V')) by taking the second 
derivative with respect to A = aeiV — V), 

{SG{V)SGiV')) = -{aef^^mV)8I{V')). 

By expressing the current in terms of Greens functions using the single-particle Breit-Wigner 
resonance conductance formula, the correlation function takes the formlll 

1 /Grr\2 ( \ 



{5G{V)5G{y')) 



f3 \ 2 J \dA 



v2 



r r duj [P^(r,r) + P_^(r,r)] 



In the absence of time-reversal symmetry, the variance changes only by the standard Dyson's 
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factor of where /3 = 1 for the orthogonal ensemble (in the presence of impurity scattering 
only) and (3 = 2 for the unitary ensemble (in the presence of a finite magnetic field or 
weak scattering by magnetic impurities that breaks time- reversal invariance). The diffusion 
propagator in the disk, Pa;(r,r), satisfies the following equation 

' -DV^ + 7 - icj] P^(r, r') = 6{r - r'). (15) 

It describes diffusion in the disk and is therefore restricted by the tunneling barrier ai z = 
and an insulating boundary at the cylinder surface p = R. The poorly conducting interface 
between the emitter and substrate is also modeled as a tunneling barrier ai z = L. These 
boundary conditions are expressed as 

^.^Uo = 0, 9,Pl=R = ^^ d,P\.=L = 0- (16) 

The correlation function is found by solving the diffusion equation, Eq. ([I5|), in the presence 
of the boundary conditions, Eq. (p^G]). The angle is set to zero without loss of generality 
and we find 

p (^j. r) = ^ '^m {cx.nmP/ R) 

nR^L (1 - mV alJ (a„^) 

iDalJR' + 7 - ^^) ' ^^^^ 
where m = 0, ±1, ±2, . . . and Jm is a Bessel function of the first kind of order m. For a 
given m, the numbers a„m are solutions of the boundary condition at the cylinder surface 
p = R, 

which may be expressed as 

We solve this boundary condition numerically in order to calculate the propagator, giving 
the variance and the correlation function for arbitrary Ld R and < p < R. 
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In the main text Fig. |^ shows the correlation function for different Lc/R and p/R = 
(main part) and for different impurity positions p/R (inset). We present here for com- 
pleteness some further numerical results. Fig. |T3| shows the calculated variance as a 
function of h'-f/T with the resonant impurity at the center of the disk p/R = and different 
values of L^/R. The short dashed line is the QOD analytic result Eq. (^ whereas the solid 
lines show asymptotics at large h'-f/T given by the Q2D analytic result Eq. (H). As expected, 
the numerical plots show behavior similar to the QOD analytic form for R (small h'y/r) 
and similar to Q2D for Lc R (large h'y/r) with a crossover at Lc ~ R. To analyze the 
effect of the position of the impurity we choose a particular value of Ly/R. The inset of 
Fig. is {5G'^) for L^/R = 1.5 and for different impurity positions. When the impurity 
position is off-center the variance has a similar qualitative form as for the impurity on the 
cylinder axis, but the fluctuations appear to be generally larger. 

Fig. |1^ shows K{AV, 0) as a function of AV/Vc and different values of Lc/R for p/R = 0.5 
(main part) and p/R = 1.0 (inset). The crossover appears to occur more slowly (over a larger 
range of L^/R) when the impurity position is off-center. 
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FIGURES 




FIG. |T]. Sketch of the resonant tunneling spectroscopy of the LDOS using an impurity 
state in a double-barrier structure. Electrons tunnel from a heavily-doped disordered emitter 
through the energetically-lowest level S of the quantum well sandwiched between the double 
barriers, so that S serves as a spectrometer of the density of states z^(-E') of the emitter. 
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FIG. H Layer structure of the double-barrier heterostructure. 
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FIG. ^. Image of LDOS fluctuations: Typical plot of the differential conductance G 
versus bias voltage ^ at i? = T and a base temperature of T = 20 mK. 
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FIG. ^. Contour plot of the differential conductance G as a function of bias voltage V 
(step 7 /iV) and magnetic field B \\ I (step 10 niT) for T = 20 niK. The excitation energy 
E on the top scale is converted from the bias voltage V, see text. 
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FIG. H The averaged differential conductance {G(B))^ of the device obtained as de- 
scribed in the text. 
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FIG. ^. Experimental correlation functions, taken at the beginning (solid, V = 10.2 mV) 
and the end (dashed, V = 13.8 mV) of the accessible voltage range. 
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FIG. ^. The correlation function at the beginning of the accessible voltage range. The 
solid line is the experimental correlation function at V = 9.8 mV and the other lines are fits 
based upon different assumptions about geometry. 
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FIG. |. The evolution of the correlation function from the beginning (bottom) to the 
end (top) of the accessible voltage range showing the experimental correlation function 
(solid) and fits for the quasi-OD (dotted) and quasi-2D model (dashed) for five voltages: 
{a)V = 13.8 mV, (b)^ = 13.0 mV, {c)V = 12.2 mV, {d)V = 11.4 mV and (e)^ = 10.6 mV. 
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FIG. ^ Theoretical form of the correlation function K{AV) as a function of AV/Vc 
with the resonant impurity at the center of the disk p/R = 0. Long dashed line is numerical 
result for Lc/ R = 1.0, dot-dashed for Lc/ R = 1.5, and short dashed for L^/R = 2.0. Solid 
and dotted lines are the analytic results in Eq. for QOD and Q2D geometry. Inset is 
K{AV, 0) as a function of AV/Vc for Lc/R = 1.5 and different impurity positions. Solid line 
is p/R = 0, dotted is p/R = 0.5, and long dashed is p/R = 1.0. 
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FIG. 0. Inelastic quasi-particle relaxation Variance of the differential conductance 

SG{EY versus excitation energy E. The inset shows the increased variance of fluctuations 
in classically high magnetic fields uJcT ~ 1. 
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FIG. 111]. Comparison of the correlation energy of LDOS fluctuations extracted from 
the amphtude [sohd hues] and the correlation function K{AV) [dashed lines] for different 
models of quasi- dimensionality. 
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FIG. Determination of the quasiparticle relaxation rate from analysis of correlation 
(dotted line) and fluctuation data (dashed line). The solid line is a fit to the theoretically 
expected inelastic particle relaxation rate 7(-E), see text for details. 
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FIG. H^. Theoretical form of the variance {SG"^) as a function of h'-f/T with the resonant 
impurity in the center of the disk p/R = 0. From the top, the long dashed line is the 
numerical result for L^/ R= 1.0, dot-dashed for L^/ R= 1.5, and short dashed for L^/R = 
2.0. The solid lines show asymptotics at large h'y/T given by the Q2D analytic result Eq. 
and the dotted line is the QOD analytic result Eq. (^). All the curves are normalized by 
the QOD analytic result Eq. (^. Inset is {6G^) as a function of h'y/T for Lr/R = 1.5 and 
different impurity positions. Solid line is p/R = 0, dotted is p/R = 0.5, and long dashed is 
p/R = 1.0. 
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FIG. |14[ Theoretical form of the correlation function K{AV,0) as a function of AV/Vc 
with the resonant impurity off center of the disk at p/R = 0.5. Long dashed line is numerical 
result for Lc/ R = 1.0, dot-dashed for Lc/ R = 1.5, and short dashed for Lc/ R = 2.0. Solid 
and dotted lines are the Q2D and QOD analytic results from Eq. (^). Inset is K{AV,0) as 
a function of AV/Vc with the resonant impurity at p/R = 1.0 with line styles the same as 
the main part. 
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